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Fluid Motions in a Sphere IV 


Thermal Instability of a Rotating Fluid Sphere Heated within 
under a Uniform Magnetic Field (2) 


By Tomikazu NAMIKAWA 


Institute of Polytechnics, Osaka City University 


(Read May 18, 1958; Received March 25, 1960) 


Abstract 


The conditions for which a fiuid sphere heated within and subject to the 
simultaneous action of a uniform magnetic field and rotation can become unstable 
via a marginal state of purely oscillatory motion are studied. Different from 
the results of the paper (Namikawa, 1957b), over-stability occurs both under astro- 


physical and terrestrial conditions. 


1. Introductions 


In earlier paper, (Namikawa, 1957c; this paper will be referred to hereafter as [1]) 
we studied the simultaneous effect of a uniform magnetic field and rotation on the 
thermal instability of a fluid sphere heated within. 

However, the discussion was limited to the case when the principle of the ex- 
change of stabilities is valid. We shall examine in this paper the possibility which 
the instability sets in as an increasing amplitude oscillation. This possibility was hnown 
to occur when rotation or a uniform magnetic field alone is present (Namikawa 1957 
a, b). 

It will be difficult to demonstrate the occurence of the over-stability or the ordinary 
cellular convection in a fluid sphere experimentally and compare the result with the 
theoretical predictions. However, this demonstration by experiment was done by Fultz 
and Nakagawa (1954-1957) in the case of thermal instability of a layer of fluid below, 
and it was in complete agreement with the theoretical predictions by Chandrasekhar 
(1952a-1956). 

Thermal Instability of a fluid sphere was studied by Chandrasekhar (1952b), Backus 
(1955), and extent of these investigation to a rotating fuid sphere was carried out by 
Takeuchi and Shimazu (1954), Chandrasekhar (1957b) Bishop (1958) and Namikawa 


(1957a). 
2. An Examination of the Principle of the Exchange of Stabilities 


We shall start from eqn. (37) of paper [1]. This derivation is given briefly in the 
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appendex. The notations are the same as in [1], but for convenience the definitions 


are listed below. 


density at the surface of sphere, 

mean density of sphere, 

coefficient of volume expansion, 
thermometric conductivty, 

magnetic permeability, 

electrical conductivity, 

1/420 

radius of sphere, 

adverse temperature gradient maintained, 
acceleration due to gravity, 

constant of gravitation, 

4/3nxpGa. 

intensity of impressed magnetic field, 

ef angular velocity of rotation, 

frequency of oscillation in case of overstability, 
ha’ 


: order of spherical harmonics, 


degree of spherical harmonics, 
2Brae 

ae * 
R.: critical value of R for the onset of a ; + Tinta 
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C=ny'a5{a5 + 20-48) g4(d+28) 7? en 
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paerif2eloss(2) +24! pe (aj +2 \2R| ) 


Ea vai {a(2 -2 R)+rl, 


5 a 
and 


a 
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From the eqn. (2), it follows that if E<0, it admits a positive real root at least. 
E<0 is therefore a sufficient condition for the principle of the exchange of stabilities 
to be valid. 

Therefore, let 


#0: (4) 


This is a necessary condition for eqn. (2) does not allow positive roots and the 
corresponding solutions lead to an exponential damping of an initial disturbance. 

We have to examine whether under these conditions the complex roots of (2) is 
able to have a positive real part. If x and y denote the real and imaginary part of 
a complex root of eqn. (2), they statisfy the equations, 


(Gx-+ Ajy*—(02' + 6 Ax +3 Bx+ Oy? +a 4+ Ant Bx’ Cx? Dx k= 07, (pu) 


and . 
yt—(10x?+4Ax+ B)y?+5at+4Ax+3Bx2+2Cx+D=0. (6) 


These are the second order algebraic eqns. for jy’; we get form (5) and (6) the 


next equations, 
Dat ada IEA DTA AB SC) AC eae sat (7s) 
and _ * ‘ ve cha - ker 
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and the eqn. (7) admits no real positive roots. Therefore, under the conditions (4) 
and (9), the eqns. (7) and (8) admit a real positive root when the inequalities (10) and (16) 
hold and by letting AD—E and AB—C->0 through negative values, we shall have the 
case when a complex root of eqn. (2), having a positive real part, tends to a purely 
imaginary limit, 

Therefore, if 


&£>0, (14) 

B>0, (15) 

AD—E<., (16) 

and AB-—C<0, (17) 


the principle of the exchange of stabilities can not be applied and we shall have a 
case of overstability. According to (3), the conditions (14), (15), (16), and (17), are 
equivalent to 


=<R, (18) 
(2+0)aj+2. Pare 72 ie (19) 


2 4 
j a; 


{21+20)43 + fo 4 AC r<{ 2 o+c420mj}2R, (20) 
5 a5 


‘D 
and 
2(14+8)a3 +2" Q<2R, (21) 
5 ag : 
The inequalities (18), (19), (20), and (21), can be combined in the form, 
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—5(20+ 3)a§+2(1-0)a2Q+T?2>0. (29) 


These four eqns. are the necessary conditions for instability to arise through os- 


cillation of increasing amplitude: it is therefore the conditions for the principle of the 
exchange of stabilities not be valid. 


If T=0 in these eqns, 
é>1 (30) 


becomes a sufficient condition for the principle of the exchange of stabilities to be valid 
and this is equal to the inequality (63) of the paper (Namikawa, 1957b). But, for T+0, 
we have no simple condition as (30) for the principle of the exchange of stabilities to 
be valid and the overstability is able to arise even if o>1. 

When the conditions (26)-(29) are satisfied, we have a case of overstability, then 
the frequency of oscillation at marginal stability will be given by from (5) and (6) 


pe Bas BAD Je (31) 
ENP 


and the critical value of R for the onset of the overstability can be obtained, letting 
AB—C+0 AD—E->0, as follows, 


4 U+ey pe 
1+26)0i+—£ ~Q aR 7 
eet fee = =«fa+or+ a}, (32) 
(14+20)1+2 529 


where a, is the first root of the eqn. Jns:/2(aj)=0. 
We can nate a relation between Q and T from (32) 


8 __ge , $8049 45(14-28)ai-T*= 0. : (33) 
5(1 aay "146 
This eqn. means that purely oscillatory motion occurs at a fixed satel of angular 
velocity of a sphere for a given intensity | of eve ee fe Ore vice fversa-civies 
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We can write the above eqn, (36) as follows ~ 


where 
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and 
A2x‘ + Box? + Cyx2+ Dox +E, =0. (40 ) 


As A, is positive and £, is negative, eqn. (39) and (40) admits a positive real root 
at least when 


E,<0. (41) 


Therefore E,<0 is a sufficient condition for the eqn. (37) has the real positive roots 
and the disturbance exponentially increases with time. There is no steady state, as FE; 
is always negative for all values of m, Q and R. 


Yate SQV (42) 


is a necessary condition for the eqn. (37) does not allow positive roots. The complex 
roots of eqn. (37) must satisfy the next two eqns. where x and y are the real and imagi- 
nary part of complex roots. 


x°+ A,x*—(10y?—4A,y—B,)x3—(6A,y?—3B,y—C,) x? 
+(5y*—4A,y3 — Cy’ + Doyt D,)x+ Ayy'—Boy'—C.y?+ Doy+k,=0 , (43 ) 
and 
(Sy— A)x*+4(A,y—Be)x?—(10y?—6.A,y?—3Bi y+ C,)x? 
—(4A,y°—3B,y’—2C, y+ D2)x+y—Apy'—Biy+C,y’+D,y—E,=0. (44 ) 
It is very difficult to examine whether the eqns. (43) and (44) are able to have real 
positive x under the conditions (42). If x=0, we have purely imaginary roots of eqns. 
(37), which satisfy the next eqns. ; 
A,(yt—A2,y®)-—CQy’+ Doyt+ £,=0. (45 ) 


and 


y—Apy—B,y+O,7+D,y—E,=0. ( 46 ) 


As the eqn. (45) is a fourth order algedraic eqn. it can be solved algebraically. 
Therefore we can get the value of R replacing in (46) with the value of y obtained 
by solving (45). But it is troublesome, so we shall return to it in another paper. We 


discuss only axially symmetric case hereafter. 


3. Examination of Onset of Instability in the Axially Symmetric Case 


From (33), we find that the overstability occurs along the line, expressed by the eqn. 
(33), in the (7, Q) plane. The critical value of R for the onset of ordieary cellular con- 


vection is found, letting E>0, (cf. paper. [1)]). 


ar-s(2o8) 2 


In order to examine the intersection of RY" and Re curves we must insert the 
relation (33) in (47), thus we have 


as [Arde tm o a0) g + 52a (48 ) 
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This has a minimum at | 
J 
= Qe 5( la Oe) eee 49 ) { 
Q=Qn=5(L+ at, [gee ( 
From the eqns. (32) and (48), we get next eqn. which is satisfied by the intersecting . 
point Q* of the RY" and Re curves. . 
{(1+6)?(1— 26?) + 69} Q¥?+.5(1+0) {6(36-+1)—2(1+6)*} hQ* + 25(1+0)2(1+2d)ai=0. (50) 
The roots of the eqn. (50) are 
[FO NORA) (51) 
where 
f,(6)=(1+6)?(1— 26?) + 68 = — 204—308—0°+20+1, 
F(0) = —6(38 +1) + 2(1+ 6)? = (20+ 1)(6?-+6 +2) >0 | 
and f3(6) = (26 +1)(0? +6 +2)? —4(— 264 — 368 — 0? — 26 +1) =0(20* + 136° + 246° + 178+ 4) >0 
(52). 
From (52), we can see that the value of @, at the intersecting point, Q*(>0), of 7 
R°" and R& curves depends on the sign of /,(0). 
Ff, has a real positive root 6,(=0.8) for a<e<d : 
i) fi>0, for 0<d<6, and R¢" and R? curves intersect at two points. 
The convection curve starts on the R,— axis at the point 
R&>c0 ~—(Q>0), = (53) 
and for Q>co becomes asymptotic to the line , 
roe ” 
2 em oe IO OY Be: aby ey, . (54. 
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Q>Q?. ! (59 ) 


This means that arising as the over-stabillity for Q@<Q*, the thermal instability 
arises as the ordinary cellular convection for Q*<Q<Q* and again arises as the over- 
stability for O@>Q*; this is what happens under astrophysical conditions (d<1). 


li) f,;=0, for d=d, and the over-stability curve, staring beneath the convection curve, 
intersects at 
5(1+¢0)ai ; 

QF a2 Ie 60 
+042 sg 
with the convection curve, lies adove the covection curve for Q>Q*, and becomes 
asymtotic to the convection curve. This means that the instsbility, arising as the over- 
stability for @<Q%, arises as the ordinary cellular convection for Q>Q* and the over- 

stability and the convection occurs at the same time for the large value of Q. 


ili) #;<0, for d>6, and the over-stability curve, starting beneath the convction curve, 
intersects at ; 


Qe OE VFR: |; ae 


and lies above the convection curve for Q>Q#. 

This means that the instability, arising as the over-stability for Q<@? and arises 
as the ordinary convection for Q>Q%. 3 

This is what happens under normal terrestrial conditions (0>1). - 


- 


iv), The two extreme conditions 60 and d>o. ; 
For 6>0, two curves become asymptotic to the lines 


Renal 437 
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and the value of Q at intersecting point tends to zero. 
In this case, the istability always as the ordinary convection. 


4. Concluding Remarks 


We examined in the paper [1] and in this paper the manner of the onset of thermal 
instability of a fluid sphere under the simultaneous action of a magnetic field and ro- 
tation. On accound of the mathematical difficulities in calculation, only the axially 
symmetric instability was studied in detail. 

In paper [1], assuming the validity of the principle of exchange of stabilities, we 
studied only the steady convection, while in this paper we examined the principle of 
exchange of stabilities and studied a marginal state of purely oscillatory motions. 

No simple sufficient condition as the inequality (63) in the paper (Namikawa, 1957b), 
6>1, can be obtained, so we have a case of over-stability both in the astrophysical 
conditions (6<1) and in the terrestrial conditions (6>1). A marginal state of purely 
oscillatory motion does not occur at arbitrary values of Q and T and does at a fixed 
value of angular velocity of a sphere for a given intensity of the applied magnetic field 
or vice versa. 

The convection curve has a minimum at the value of Q given by (49), while the 
over-stability curve is a monotoneous increasing with Q. In the case of the paper 
(Namikawa, 1957b), the over-stability curve always starts above the convection curve, 
but in the case of this paper the convection curve always starts above the over-stability 
curve. This meavs that the over-stability occurs at the small value of Q. 

a) When 7 <0=0.8, the over-stbility occurs at the small (Q<Q#) or at the large 


(Q>Q#) values of Q, while the ordinary convection arises for the moderate values of 
Q@ given by (58). Under astrophysical condition, for example, in the so called hydrogen- 
convection zone in the solar atmosphere (Chandrasekhar, 1952a) e~10'%cm2/sec, 7~7X 
10®cm?/sec and 6~7X1077 <1. : 
This is the case above mentioed. _ = 
b) When % pare ro the over-stability occurs for the values of a given by (ae or (61) 
and for the ice, values of Q, the ordinary convection occurs. 
_ Under terrestrial conditions, for ex , 


* 
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field. 
Recently, Rikitake (1959 a, b) studied thermo- -magneto-hydrodynamic oscillation in 
the Earth’s core under the steady Poroidal and Toroidal magnetic field. He got twenty 


years for the period of the oscillation, which is in accordance with that of a part of 
the geomagnetic secular variations. 


5. Appendix 


1. Equation of motion 

A brief derivation of the eqn. (1) is given below. 

Consider a homogeneous electrical conducting fluid sphere of radius a rotating with 
an angular velocity about Z-axis under a uniform magnetic field, the direction of which 
is parallcl to Z-axis. A heat generating source is supposed. The rate of heat generation 
is assumed such that in the absence of conduction and convection the temperature 
would rise at a uniforme rate «. The equation of motion, of Maxwell, and of heat con- 
duction appropriate to the problem on hand are 


p Ou 


or —=+o(u-p)u=—pp+2oux 2+ opus op V+ nx, (AwT) 
90 4 div(ou)=0 A. 2) 
a +div(ou)=0, ( 
curl H=47j, (A. 3) 
cade = 2, 2 (A. 4) 
at 
div H=0, - (Ase) 
> div E=4rzc’q, (A. 6) 
j=oF+uxH), (A. 7) 


= ss +(u-7)T=Kp°T+e, (A. 8) 
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where 0 and A are perturbation in temperature and magnetic field. 
Taking the curl of eqns. (A.10) and (A.11), we get 


Om Ou wH oO 
6xr+29-24 lh+vp2 , (A.13) 
“Ot sa ulapinil 0z =F Ap) Oz Pralieoua exkaie 
and 
© curlhka=Ho® +npcurlh). (A.14) 
ot Oz 


Again taking the curl of eqn. (A.13), we have 


HT 0 pep yyy, (A.15) 


0 (pe 
ee = l(voxr)— 
rnb u)=ycurl(vd xr) tno” 02 


where @=curlv. 


Multiplying eqns. (A.11) (A.13) (A.14) and (A.15) by r, we get 


2 h,)= =Hroh +07 X(rhy) , (A.16) 
LF (r0,)= 2.0r Ou OH! m ae MN 2 (curl) +yr(r0,) (A.17) 
ates +yp?r(curlh), , (A.18) 
Oz 
and 
UP ese Pe 20r 9@ Balt Bs Hy é ‘bolgeacin wees (A.19) 
ot 0z Arno car 
where 
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As in [1], we shall examine the first type motions taking only the first four terms 


in consideration. 
There is no second type motion because the continuity of the magnetic field can 
not be satisfied on boundary surface. 
m m+. ‘mm m 
Un > Vn Tm Smet and@ ” by the spherical Bessel 
r Yr 


, , 
r2 yr? 


We expand five scalars 


functions | el where a,’s is the roots of Jns:/2(a;)=0. With above expansions 
r 
the boundary conditions are satisfied. 


If all quantities vary as e*’, substituting the five scalars, expanded by the spherical 
Bessel functions, in the basic eqns., we get the eqn. (1). 


6. Acknowledgments 


The author wishes to express his thanks to Prof. S. Chandrasekhar for his discussion 
of my previous papers and for informing me of Bishop’s paper (1958). My thanks 
are due to the Director M. Hasegawa of Fukui University and Prof. Y. Saito of Osaka 
city University for their continuous interest and encouragements in the course of this 
study and to Dr. T. Rikitake of Tokyo University for his valuable discussions and advice 


on this problem. 


References 


Backus F.E. (1955) Phil. Mag. 46, 1310. 

Bishop F.E. (1958) Phil. Mag. 8, 3, 1342. 

Bullard E. (1948) M.N.R.A.S. Geophys. Suppl. 5, 248. 
Bullard E. and Gellman H. (1954) Phil. Trans. 247, 27. 
Chandrasekhar S. (1952a) Phil. Mag. 7 48, 501. 
Chandrasekhar S. (1952b) Phil. Mag. 7, 43. 1317. 


~ 


Chandrasekhar S. (1953a) Proc. Roy. Soc. A 216, 293. a 
racsekhat eS ten sea a Ae ean i: Re ts) Sih) giz Ls bie = avi sa FS ° 


Magnetic Viscosity of Magnetite* 


By Yoshio SHImizu 


Geophysical Institute, Tokyo University 
(Read October 17, 1959; Received April 14, 1960) 


Abstract 


Magnetic viscosity of ensemble of medium-sized grains of natural magnetites 
was examined at various temperatures. It was found that magnetic viscosity 
coefficient S, defined as J—J,=S(Q+logt) where J and J, are intensity of magne- 
tization at time ¢ and ¢=0, is proportional to external magnetic field, and S isa 
linear function of temperature, except at temperatures close to the Curie point 
and — 160°C, so far as the Rayleigh region of magnetization is concerned. It 
seems also that S tends to a finite value according as grain size tends to the order 


of a single domain. 
For ensemble of magnetite grains, S in expression of 4J=S(Q+logt) is very 
small, and the limit of half life time c for change in thermoremanent magneti- 


i, 
ears. 
2S VA 


Further, the half life time for sedimentary rocks, in which directions of grain 


zation, expressed by t= 


magnetization are scattered, with reduction factor w, is about 10!° years for w=1/5. 
These results may show that remanent magnetization of igneous and sedimentary 
rocks have been stable against thermal viscosity during a geologically long time. 


1. Introduction 


It has long been ‘known that some igneous rocks or baked sediments have very 


stable remanent magnetization, which was fixed in the direction of the earth’s field is 
earns) S ; : : 
during the process of cooling from above the Curie temperature (Konigsberger, 1932 ; i= 


Thellier, 1938; Nagata, 1940). Since this permanent magnetization was discovered, we 
ve expected that the thermoremanent magnetization of rocks may have fossieds 
direction ei tee gree pat at the age of et ttt has” also 


netization Ww 
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the magnetic field before folding (1949). On the other hand, Kawai observed in certain 
sedimentary rocks that remanent magnetizations of such rocks are always shifting 
their directions towards the ambient magnetic field when the rocks were replaced from 
one direction to another in laboratory (1953). More recently occurrences of such unstable 
remanent magnetizations have been reported (Creer 1957; Akimoto and Kushiro, 1960). 
Prior to these, Thellier indicated that the intensity of the isothermal remanent magne- 
tization of some clay decreases proportionally to the logarithm of time for the time 
intervals from 20 sec to 5x10®sec. (1938). Such change of magnetization with time has 
long been known as “ magnetic viscosity” or “magnetic after-effect” in the case of 
metallic ferromagnetics, (Ewing, 1885; Rayleigh, 1887). 

It is intended, in the present study, to make investigation of magnetic viscosity of 
magnetite as the most common ferromagnetic minerals in rocks, by considering its 
relevance to palaeomagnetism. Measurements were made for various magnetic proper- 
ties, i.e. susceptibilities, coercive force, magnetic viscosity and thermoremanent magne- 
tization and their thermal variations within the wide range (90°~850°K) of tempera- 
tures. The results of the measurements will be given in § 4 and in the last section 
possible mechanism of the observed magnetic viscosity of remanent magnetization will 
be discussed. 


2. Specimens 


Specimens used in the present investigations were natural magnetite grains of iron 
ore from Kamaishi Mine, Japan. It was confirmed by microscopic examination that 
each grain im situ was composed of crystallites. The chemical composition, lattice 
parameter, and intrinsic magnetization of the specimens are summarized in Table I. 

Table I. Chemical composition, lattice parameter and spontaneous 
magnetizatoin of specimens. : 


| weight per cent mol per cent 


— iH so Sie i ——— = SS ee Fe = = 
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Table II. Summary 
present investigations. 


of structure sensitive magnetic properties of magnetite grains in the 


Ree acne es differential susceptibilities : 
SS | - “size | reversible. ] irreversible f pees 
| a aR “| | in emu/gr | in emu/gr mortaal we 
1 270 0.698 | Giogites Pt SOiyama 1 
2 100 0.493 0.190 19 
3 ~2 0.429 0.132 50 
4 ~1 0.314 0.125 66 
viscosity viscosity thermoremanent 
_ coefficient field constant magnetization 
in emu/gr x 10-2 in Oe in emu/gr 
0.85 0.096 5,54 
0.66 0.035 3262 
0.60 0.046 Goo 
0.49 0.039 — 


3. Method and Apparatus 


The specimen was mounted in a non-inductively wound and electrically heated 


tublar oven for heating experiments 


up to 600°C. In cooling experiments, a Dewar 


vessel containing either liquid oxygen or “dry ice” moistened with alcohol was used. 


External magnetic field was applied by means of a water-cooled solenoid which is 30cm 
long and 5cm in diameter. The required stability of the magnetic field was obtained 


by supplying the current to the solenoid from a 6v. battery. 
Measurements of the intensity of magnetization were made by an astatic magneto- 


meter. Arrangements of the apparatus are schematically shown in Fig. 1. 


Fig. 1. Schematic view of arrangement of ap- 
paratus used for measurements of magnetic 
susceptibility and viscosity. In the figure, the 
system SMN denotes an astatic magnetometer, 
F and F’ are field coils, J is a specimen, C is 
a compensation coil, and T is a thermo-couple. 
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Temperature meaurements were made by a calibrated thermo-couple having one 
junction placed near the mid-point of the specimen. It was proved that the temperature 
of the specimen could be maintained for several hours at any constant temperature 
between 15°C and 400°C within 2°C in accuracy and a little worse at higher tempera- 
tures; with solid carbon dioxide in the Dewar vessel the temperature was —77°C+1°C 
for several hours, and when liquid oxygen was used the temperature was maintained at 
—183°C-41°C for the same length of time. For the purpose of giving intermediate 
temperatures between —77°C and —183°C, the method illustrated in Fig. 2 was used. 
By this method, however, it was difficult to maintain the constant temperature for more 


Fig. 2. Schematic figure of 
apparatus for temperature 
control between 90°K and 
196°K. 


than several minutes. 

For measurements of magnetic viscosity for a point on 
the virgin curve of magnetization, the following procedure 
was adopted. The specimen was carefully demagnetized 
thermally in a non-magnetic space and was brought into the 
apparatus, in which the vertical component of the terrestrial 
magnetic field was kept compensated by appropriate current 
of the field solenoid. Then the whole apparatus was allowed 
to reach thermal equilibrium. The deflection of the magne- 
tometer corresponding to instantaneous magnetization due 
to the application of the external field, was compensated by 
appropriate currents through the carlibration coil, since we are 
interested only in the subsequent small variation of magneti 
zation with time. During the measurements, no change of 
the currents in the field coil and in the calibration coil was 
detected. 

The measurements of the total, reversible and irreversible 
differential susceptibilities were made, by observing small 
changes in magnetization caused by slight changes in applied 
field. The definitions of these quantities are 


cD Le re 
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Fgi. 3. [Illustration of measurement of reversible and 
irreversible differential susceptibilities, where «Hs denotes 
a small variation of magnetic field. (s=0,1 and 2). 
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applied magnetic field H in the present experiments is 3.3 Oe. The result shows that 
apparent time variation of magnetization of magnetite grains, 4/, is proportional to 
logt within the time intervals from 30 to 2,000 seconds, that is, 


Al=const.+ Slogt, (ed. 5) 


where the slope, S, is called the magnetic viscosity coefficient. If the applied field 
H (=3.30e) is suppressed after an instant application, the residual magnetization 4/, 
decreases gradually. (4/,, logt) curve shown in Fig. 5 represents 4/, as a linearly 
decreasing function of logt, the slope of which is about half of that of the above- 
mentioned experiment; i.e. S of the former is 0.45x10-® emu/gr and that of the latter 
is 0.85107? emu/gr for the specimen No. 1. 


$= 045 * 107 
15 


Specimen No.1 


AJ (emu/grx 107) 


{00 om it 10°. 
Time t(sec) 


Fig. 5. Variation of remanent magnetization as a logarithmic pee of time ¢ after a 
_ suppression in ahaa Bese field. — = ' 
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Applied Mognetic Field Coercive Force 
Fig. 6. Variation of magnetic viscosity Fig. 7. Variation of magnetic viscosity coefficient 
coefficient S as a function of the applied S as a function of reciprocal value of coercive 
magnetic field intensity. force. 


Fig. 8 illustrates the temperature dependence of the magnetic viscosity coefficient 
S as deduced from the experiments shown in Fig. 4. The applied field is 3.3 Oe through- 
out all three experiments. Measurements near the teinperatures of liquid oxygen were 
done only for the specimen No. 1. At the temperatures near the Curie point and near 
—160°C, the experimets required extremly highly accurate temperature control. Even 
under a good control of temperature, the duration of observation had to be shortened 
to only several minutes. In these figures, the most remarkable characteristics are 

(i) There are two high peaks of S at temperatures near the Curie point and near the 
transition point at —160°C. These peaks will be refered to as S, and S;, in the following. 

(ii) Between the two peaks S, and S, within a fairly wide range of temperature, 
S varies as a linear function of absolute temperature 7. 

(iii) The height and the temperature of the peak S, depend on the grain size of 
the specimen, and hence on the coercive force. 

In order to determine the temperature tendency of viscosity field constant Sy=S: X;,,, 
which has been introduced in the Néel’s theory given below (1950, 1955), the temperature 
variations of reversible and irreversible differential susceptibilities are measured between 
196°K and 850°K. An example of the results is shown in Fig. 9. 

Fig. 10 (a), (b) and (c) show the thermal dependence of S, as deduced from the 
experiments in Fig. 8 and 9, together with the temperature dependence of the coercive 
force. The figures indicate that S, varies linearly within only the limited range of tem- 
peratures from 200°K to 400°K or more, and at about 600°K it reaches the broad maxi- 
mum, and then shows a minimum. Above this minimum, the S,—T curve has a sharp 
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apex near the Curie point. Simi- 
larly as in Fig. 8, the higher the 
coercive force of specimen is, the 
nearer is the apex to the Curie 
‘point. Such behavours of Sy will 
be discussed in next section. 

In Fig. 11 are plotted the 
values of total, reversible and 
irreversible differential suscepti- 
bilities against the intensities of 
the applied magnetic field. The 
figure shows that the reversible 
susceptibility is independent of 
the field H, and the irreversible 
one is proportional to H. Hence, 
the instantaneous magnetization 
after the application of the 
magnetic field H may be ex- 
pressed empirically as follows 


(5 ) 

where A denotes the reversible 
susceptibility and 2BH is equ- 
ivalent to irreversible differenti- 
al susceptibility. Eq. (5)  ex- 
presses essentially the Rayleigh’s 
first law of magnetization 
weak field 
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° ry 0.4 3 
i ee Fig. 11. Variation of total, reversible 


. “hg ‘ a - and irreversible differential suscepti- 
bilities as a function of applied mag- 
i netic field. 


given by 
A,(t)=Sy(Q+logt) , (6.9 


where Q is a numerical constant of the order of 40 or 50, ¢ the time interval since the 
beginning of the application of the field H, and S,, which has dimension of magnetic 
field, is a characteristic constant of the specimen, depending on temperature. Thus, 
intensity of magnetization J after a time ¢ has been derived by Néel as 


T=1h+%irrSv(Q+ log?) , ( 7 ) 


where J, is magnetization in case of absence of the viscosity field and y;,, the irre- 
versible differential susceptibility. 

In the case of Rayleigh region, the instant magnetization J, is equal to AH+ BH* 
and the irreversible susceptibility y,,, is equal to 2BH as can be seen in $4. Thus, 
eq. (7) can be expressed as 


I= AH+ BH*?+2BHS,(Q+log t). (8) 


3 After suppression of the applied field H, the remanent magnetization is given by 
+BH 2 and the irreversible susceptibility by —BH, so that the remanent magnetization 
I, at the time ¢ after the suppression of the field becomes 
sh * 


. ; “Lei oomiseal 


1, = 7 BH? — BHS(logt—logt’) , bias: ces (9) 
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with decreas eee 
ease in grain size but the former tends to a limit value for a single domain 
value as its extremum condition. 
Now, quantity S,, which is j i 
; q é y Sy, which is introduced in eq. (6) is significant in relation to mechanism 
Or magnetic viscosity. Accordine ee ; 
gnetic viscosity. According to the Néel’s interpretation, the Bloch wall may cross 
obst s by f the inter ispersion fi ich i 
acles by help of the internal dispersion field which is derived from small oscillations 
9) o - me Fe ake 5 Ai 
of each spontaneous magnetization around the mean direction of domain magnetization 


(1950, 1951 and 1955). Based on the above-mentioned model, it is given that 


Se AnkT Qu1/2 
=,/® ee (10) 


where & is Boltzmann constant, T absolute temperature, and v is volume of Barkhausen 
discontinuity. 

On the other hand, another theory has been proposed by Street and Woolley (1949, 
1950, 1952 and 1955). In their theory, the thermal agitations are supposed to affect 
heights of the barriers which opposed the wall movements. Using the theory of rate 
process, the following relation was derived, i.e., 


kT 
= Qui,” ee 
where /, is spontaneous magnetization and v is volume of the Barkhausen domain. The 
Néel’s theory on magnetic viscosity for single domain grains is essentially similar to” 
the Street and Woolley’s (Néel; 1949 and 1955). 

As can be seen in eq. (10) and eq. (11), the thermal variation of S, is given differ- 
ently in the two theories; if v is independent of 7, S, varies proportionally to T'/? in eq. 
(10), but to T in eq. (11). 

In the results of the present experiments, however, Sy showed some complicated 
features, as can be seen in Fig. 10 (a), (b) and (c). However, by refering to the above- 
mentioned two theories, the present results may be interpreted for three temperature 
ranges as in the followings. 

(i) ca 200°K<T<ca 400°K. In this range of temperature, Sy is proportional to 7. 
According to the theories summarized above, the viscosity field in this range is not 
likely to be controlled by the internal dispersion field, but by the fluctuation of the 
heights of barriers. 

(ii) ca 400°K<T<ca 700°K. Sy increases slowly, deviating from the linearity of 
S,-T relation, and reaches maximum and then decreases with increase in temperature. 
On the other hand, S=7;,,-Sy is about a linear increase function of T, as can be seen 
in Fig. 9, whereas 7%; augments gradually with temperature. Hence, the increase of 
S, with temperature becomes slow and deviates from the linear relation between Sy 
and 7. 

(iii) ca 700°K<7. At temperatures very close to the Curie point, with increase 
in temperature S, increases again rapidly, forms a very sharp apex and then tends to 
zero rapidly. The temperature at the apex of Sy coincides with that of S and also that 


of the maximum of 7. This sharp inarease of S, may prove that the thermal increase 
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of S is now more rapid than that of y;,,.. Further, the temperature of the apex is differ- 
ent from specimen to specimen, being dependent on grain size. It is also concluded that 
coercive force of the specimen becomes, at this temperature range, as small as the order 
of magnitude of external magnetic field. This result might be interpreted as follows. 
At the temperatures very close to the Curie point, crystal anisotropy of specimen becomes 
very small and thickness of Bloch wall becomes larger and larger until the meaning 
of domain conception is lost, whereas the shape anisotropy does not decrease so rapidly 
compared with the crystal anisotropy, because the shape anisotropy energy is proportion- 
al to J? and crystal energy, according to Zener-Carr’s theory, is proportional to the higher 
power than J? (Zener, 1954: Carr, 1958). In such a state, each grain may be regarded 
as a single domain grain which has a critical field H,=K/I, for rotation of spontane- 
ous magnetization, where K denotes shape energy. Magnetic viscosity of such single 
domains may be treated by the Néel’s theory for single domain grains. That is to say, 
provided that intensity of magnetic field H is close to critical field H, of grains, S and 
Sy, are approximately given by 


2kT 
Se . (12) 
pes Ee 8) 
Hk l: (13) 


i Opet Vien a) : 


where H,,, denotes maximum critical field, and var is maximum volume of single domain 
grains. If H,m approaches to H with increase in temperature, factor (H.,—H) or Hem/ 
(Hom—H) in eq. (12) or eq. (13) becomes rapidly large. 

The present results are compared with those of metals and alloys made by Courvoisier 
(for steel ; 1954), Barbier (for alnico ; 1953), Street, Woolley and their colleagues (for alnico 
and other several alloys; 1949, 1951 and 1955), and Pescetti and Barbier (for alloy of 
Cu, Ni and Be; 1956). Generally speaking, the present results are in qualitative agree- 
ment with those for metallic ferromagnetics, though the apex of S, at a temperature 
near the Curie point is not found in their results. 


6. Stability of Remanent Magnetization 


Viewed from the results obtained in the present study, stability of remanent magnet- 
izations obtained by various procedures of magnetization such as thermoremanent 
magnetization (TRM) and detrital magnetization (DRM), etc. will be discussed. In 
ferromagnetic specimen having such permanent magnetization, the Bloch wall has been 
trapped by an obstacle at a position different from that corresponding to zero magnet- 
ization. Under the constant temperature and constant pressure, distribution of such 
obstacles, which are due to imperfection of crystal such as fluctuations of internal stress, 
impurities, or shape anisotropy etc., is expected to be kept constant in the interior of 
the crystal, provided that there is no secondary chemical alteration. Under the effect 
of thermal agitations, the Bloch wall may cross over trapping obstacle. Even if this 
probability is negligibly small in the laboratory time scale, its effect should become 
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appreciably large for very long time such as geological time scale. For the laboratory 
time scale, it is proved that the relation dJcc log ¢ holds for natural magnetites also. 
However, it may be clear, on the other hand, that the total amount of change in in- 
tensity of remanent magnetization caused by magnetic viscosity should tend to a finite 
value when ¢ becomes infinite. Therefore, amount of decay of remanent magnetization 
estimated by means of the law JJ cc logt during a certain time interval should be upper- 
most limit, and consequently the decay time for a certain amount of change in /, 
calculated by the formula should be its lowest limit. By considering the above-mentioned 
condition, stability of remanent magnetization for geologically long time could be 
estimated as in the following. 

When the direction of the external magnetic field is antiparallel to that of ther- 
moremanent magnetization of a rock specimen, change in magnetization Xt) during ¢ 
is given by 


[(t)=I,,—S(logt—logt’) , (14) 


where /,, is the initial thermoremanent magnetization and ?’ is the time intervals when 
the magnetic field kept parallel to the direction of thermoremanent magnetization at 
room temperature. If ¢ is negligibly small compared with ¢, the “half-life” time of 
thermoremanent magnetization will be given by 
ly 
MO or (15%) 

Numerical estimation based upon actual values given by figures in Table II, results 
in that the “half-life” time of the thermoremanent magnetization of ensemble of mag- 
netite grains is in order of 10!2° years, which is much larger than the age of the earth. 
Thus, it may be concluded that thermoremanent magnetization of igneous rocks is 
extremely stable against thermal viscosity. In other words, thermoremanent magnet- 
ization of igneous rocks ought to have well reserved its initial condition during their 
whole life, provided that there has been no appreciable chemical modification. 

In the case of detrital remanent magnetization of sedimentary rocks, which is con- 
sidered to be due to statistical alignment of magnetite along the earth’s magnetic field 
whose remanent magnetization is originally due to thermoremanent or chemical rema- 
nent magnetization (Kobayashi, 1956). The “half-life” time may be given as 


Wir —Slogt, (16 ) 


where w denotes an efficienency of the alignment. 
Eq. (16) gives, for example, that ¢~10* years for w=1/3 and t~3 x 105 years for w= 


1/10. According to the results of evaluation of w of some actual sediments, its magni- 
tude is between i and + (Nagata, 1953). In such a case, stability of the detrital 
remanent magnetization against thermal viscosity is fairly high. 

On the bases of above mentioned estimations, we can place some reliance upon the 
results of palaeomagnetic studies by means of igneous or sedimentary rocks, provided 
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that due cautions are taken in selecting the rocks. 


rh 
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The Structure of the Atmospheric Electric Field 


By Toshio OGAWA 


Geophysical Institute, Kyoto University 
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Abstract 


Simultaneous records of the atmospheric electric potential gradient at three 
stations which make a regular triangle with sides of about two km, were compared 
with each other during the period of January 8-13, 1956, including a cloudy day 
as well as fine days. Expanding the daily course of the potential gradient on 
each day and at each station into Fourier series of 12 terms, the spectral distri- 
bution of the square of the amplitude of each component was obtained and the 
spectrum of the field was given by m~‘ in an undisturbed day, m expressing the 
harmonic term. On theoretical grounds, assuming the distribution of space charge 
concentration characterized by its scale and its density within the lower layer of 
the atmosphere, the spectral distribution of field energy was suggested to be 
represented by u~7 in the intermediate range of m. The result was proved to be 
consistent with the above analytical result. The factor causing the disturbance 
in the spectrum was looked for into the atmosphere using the aerological data. 
It was found that the turbulent situation of the air in the lower altitudes up to 


three km would control the structure of the atmospheric electric field. 


1. Introduction 


Horizontal distribution and time variation of the atmospheric electric potential 
gradient have been investigated in a series of surface measurements carried out in 
Kyoto since 1955. The purpose of this paper is to find out the structure of the atmos- 
pheric electric field in electrically undisturbed day and the cause of the perturbation 
in the field. 

On land areas of the earth, especially less removed from sources of pollution, the 
atmospheric electric field which is usually measured within the turbulent and contami- 
nated exchange layer is subject to great local variations, which mask almost entirely 
the universal variation of the field. Even in a so called undisturbed day there are so 
many kinds of variations, regular and irregular, in the field at the earth’s surface that 
it is impossible to attribute each variation to a corresponding origin. To see the gener- 
al feature of the field, the statistical analysis of large amount of data is usually re- 
quired, which offered much of the existing information, leaving the choice of periods 
for the analysis to the investigator. In the statistical procedure of analysis, irregular 
ular variations are emphasized, and the charac- 
change with local time is obtained, which is 
pheric electric field. However, it is not 


variations are eliminated and only reg 
terizing and pronouncing daily curve of 
regarded as the representative of the atmos 

(139 ) 
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enough to conclude the representativeness of the variation curve of the field thus ob- 
tained to that in fair weather. Can it be acceptable that the eliminated irregular 
variations with periods smaller than diurnal and semi-diurnal are not significant in the 
atmospheric electric field profile? Such irregular variations must be related with the 
turbulence in the atmosphere, and they are the substantial components which the field 
consists of. 

As has been pointed out by Kasemir (1950), atmospheric electric processes such as 
changes in the potential gradient and space charge density are related with the relaxation 
time of the air. So far as treating variations with rather larger periods than the 
relaxation time, the atmosphere is considered to be in a state of quasistatic equilibrium 
and the principle of the quasistatic state can be adopted. Since actual values of the 
atmospheric relaxation time are from 5 to 40 minutes at the earth’s surface, it is possi- 
ble to apply the quasistatic principle to variations with periods of over one hour. 

In this state the vertical current density 7 is considered to be constant with altitude, 
and the space charge density o which is arised from the conductivity change with 
altitude is given by 


p=eiZ(_), (1) 


where « is the atmospheric permittivity and 2 is the conductivity. As the conduc- 
tivity of the atmosphere increases almost exponentially with height, the space charge 
density is considered to decrease with height in the same manner. In actual, smaller 
conductivities due to the condensation nuclei concentration produce the positive space 
charge in the lowest altitude of the atmosphere. 

The intricate variations in the atmospheric electric field at the earth’s surface are, 
in fair weather, considered to be caused by changes of the space charge configuration, 
which is, of course, influenced by the atmospheric turbulence. It is also the purpose 
of this paper to learn how the atmospheric electric field is modified by the day-to-day 
variation in the intensity of the turbulence in the lower atmosphere. 


2. Experimental and Analytical Results 


The circumstance of the measurement for the atmospheric electricity is as follows. 
Measurements were made at three points which make a regular triangle with sides 
of about two km in the north-east district in Kyoto city which has the population of 
about one million and no big industrial area. The reason why the three points were 
used is to find out the chance that very local changes in specified points might disturb 
the generality of the results. The surface stations were all located atop the buildings 
nearly 10-15 m high to avoid the disturbances near the ground. Throughout the year 
and especially in the cold season, it is excelled in westerly wind over Kyoto. 

The reproduction of the records for the period of January 8-13, 1956, is presented 
in Fig. 1, the values of the potential gradient being normalized to the mean value for 
each day and each station. According to the report of Kyoto Meteorological Observa- 
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POTENTIAL GRADIENT 


fe) l2 24 12 24 12 24 12 24 12 24 12 24 
LOCAL TIME 

Fig. 1. Variations of the potential gradient at Yoshida (upper), Sokokuji 

(middle) and Shimogamo (lower) on January 8-13, 1956. The value of the 

potential gradient is normalized to its mean value for each day and each station. 


Table 1. Meteorological conditions during the period of January 8-13, 1956 


ad 5 Date ‘ | | 


rhe Jan. 8 Jan. 9 Jan. 10 Jan. 11 Jane l2 Jan. 13 


Meteorological 
conditions 


Cloudiness, 2.0 1.0 1.8 0.0 0.8 5.3 
average, 0-10 


Wind speed, 21 2.0 1.4 1.4 0.9 0.3 
average (m/s) 


Sunlit time 6.9 8.8 ‘Say 5.1 7.6 0.7 
(hour) 


tory, since January 8, Japan Islands were in the typical winter type distribution of 
atmospheric pressure, high in the west and low in the east. The surface meteorological 
conditions—cloudiness, wind speed and sunlit time—on each day are -given in Table 1. 
During this period there was no precipitation except January 11 when a minor shower 
cloud with the precipitation 0.0 mm swept over the stations at about 17 o’clock. The 
effect of the minor shower was quite insignificant as shown with stars in the figure. 
At no time during this period was a negative value of the electric field observed. Nor 
did any other changes in the field occur which could be attributed to visible meteoro- 
logical effects. The remarkable feature given in Fig. 1 is that double periodic curves 
with one maximum in the forenoon and in the afternoon respectively were recorded 
on January 8, 9, 10 and 12 at all the stations, but on January 11 and 13 single periodic 


curves with the maximum at about noon were found. A same example was given by 
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Israél (1955), and it was examined from the synoptical aspect by Dolezalek (1957). In 
a comparison between three stations the general features with larger periods such as 
diurnal or semi-diurnal are similar to each other, but there can be seen some different 
points in detail. 

Using the hourly values normalized to the mean value, over a 24 hours period for 
each day and each station, the variations were expanded into Fourier series of 12 terms 
and the daily spectral distributions of the square of the amplitude, £7, of each com- 
ponent were obtained. The result of this procedure, which is rather divergent, is shown 
with logarithmic scales in Fig. 2 but it can be seen that there is to some degree the 
regularity in the spectra. In general, there is a peak at the second term in the spectral 
distribution and the values of £2 decrease when terms become higher, and in the range 
higher than 6th or 7th term the spectral distribution is subject to a great fluctuation. 
In the intermediate range from 2nd term to 6th or 7th term, however, rather fixed 
value of characteristic slope is seen on January 10,12 and 9. The most typical example 
(January 10) gives the spectral distribution of m-*. On the other hand the estimation 
of the value of the power is not possible on January 8, 11 and 13. The question arises 
as to what value of the power to be in an undisturbed condition,or what is responsible 
for the disturbance of the spectrum. A theoretical consideration will suggest in the 
next section that —4 is available to the value of the power of the spectrum with an 
appropriate model of the field structure when there is not any disturbance. The pre- 
ceding experimental value on January 10 is in agreement with this value. 

The factor which controls the atmospheric electric field profile is the distribution 
of pollution particles with altitude, which is determined by the intensity of turbulence 
in the lower atmosphere. Then it will be the turbulent situation in the lower atmos- 
phere that disturbs the spectrum of the field. Air temperature and wind volocity distri- 
bution with altitude may be the main determining factors in the turbulence of the air. 
The air temperature and the wind velocity at Yonago Meteorological Observatory were 
shown with altitude in Fig. 3. The values at 12 o’clock in the mean local time were 
used for the elements. Yonago is located about 200km to the west of Kyoto and the 
upper air mass trajectory for 850 mb level passing over Kyoto came from over Yonago 
during the period. 

Sagalyn (1958) gave the criterion to see the exchange layer. ‘“ The temperature 
distributions always show an adiabatic temperature gradient throughout the exchange 
layer with a change to a more stable lapse rate at the top of this region.” This cri- 
terion suggests that the upper boundary of the exchange layer is 2-3km above the 
ground in this period. The space charge constituting the atmospheric electric field 
can not leak out from the top of this exchange layer and may be distributed within 
it because of the stable region over the exchange layer. Thus the turbulent situation 
of the air was examined up to 3km. To find out the effect of the turbulence of the 
air on the atmospheric electric field, the author calculated the Richardson number 
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which represents a criterion to decide whether turbulence will subside or increase, 
where T is the absolute temperature, v is the mean velocity, g is the acceleration due 
to gravity and J’ is the adiabatic lapse rate. This is given in Fig. 3 in which the 
calculated values are taken in steps every 500m. From a comparison between Fig. 2 
and Fig. 3 it can be seen that there is a correspondency from day to day between 
the spectrum of £%, and the distribution of R; with altitude. The thicker the layer of 
R;>1 in Fig. 3 is, the more systematic the spectral distribution of £? in Fig. 2 becomes, 
and the value of the power of the spectrum tends to converge a definite one. Thus 


R= ’ (2) 
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the instability of the air may result in the modification of the spectrum. 


3. Theoretical Consideration 


As in the previous analysis, the time variation of the atmospheric electric potential 
gradient at the earth’s surface may be expressed in terms of the Fourier series: 


E=Eo+ DE, sin(nt+ap) , (3) 


that is, it is composed of various simple harmonic components superposed on the mean 
value. 

It is interesting to give one description of the atmospheric electric field by intro- 
ducing following the two concepts ; geometrical scale and electrostatic energy. In the first 
place, space charge configuration causing variations in the electric field should bec on- 
sidered. Suppose that there are many ranks of space charge concentrations by which 
the corresponding harmonic terms of the field variation can be caused, and any rank 
of space charge concentration is in a limited region with a characteristic length A 
(scale) and also electrically characterized by its density o. The scale A is proportional 
to the reciprocal of the rank m,°1.e. 


too (4) 3 


It is assumed that there is no mutual interaction between concentrations both dyna- 
mically and electrically, and each of concentrations is located as shown in Fig. 4. It 
follows that the height of each rank of concentrations, H, is proportional to the scale 


ee 


as 


Hed. r (5) 


Next, any rank of space charge concentration is assumed to have its own density 
which is proportional to its scale, i.e. 2 S - 


< 4 rd 
™ 
= .s 


ocd. . (6) 
The space oe density of the concentration with a a scale contributes to the total 
space charge density distribution with altitude in the higher levels of altitude as well — 
| as in the lower levels. On the other hand, the space charge density sritetbeitaica 
<< ibutes in only pele veli teas bey easton ae 2 
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Fig. 4. Schematic diagram showing the inter- 
relation between the scale A and the density 
op of the space charge concentration. For 
example, 2nd rank of space charge con- 
centration has the scale 4, and the density 
02. H is the height of the concentration. 
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Table 2. Interrelation between the values of , A, H and p. 
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where a region of a sphere is available for three dimentional model. Then the electro- 
static energy is given by 
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Ea 28 (12) 


Substituting into (8) and referring to (5), 


UC a (13) 
From (10) and (13) 

E*cc A, (14) 
or 

Bron . (15) 


This is in agreement with the above analytical result of January 10, 1956. 

To see the order of the value of the actual scale, 2 or 3km which corresponds to 
the height of the upper boundary of the exchange layer should be used as the value of 
the largest scale A, corresponding to the first rank. For example, as a scale of 
6th rank it results in about 400-500 m. It is seen in Fig. 3 that in the lowest 500m 
of the atmosphere there is an layer where R;<1, during the period of January 8-13 
except January 12, and the space charge configuration in this lowest unstable layer 
may be the cause of the great fluctuations at higher terms in the spectra in Fig. 2; 
that is, the space charge concentrations with smaller scales than 500m are apt to be 
influenced by the turbulent situation of the air in the lowest levels of the atmosphere. 


4. Concluding Remarks 


The analytical results discussed above have been limited to the measurements ob- 
tained during only successive six days and may be not conclusive, some matters remaining 
to be confirmed in future. However, the results deduced in this paper may have an 
important meaning for the description of the atmospheric electric field; that is, the 


atmospheric electric field will be sufficiently understood by considering its perturbation 


which may be missed in the statistical method. Therefore, just as the daily curve of 


siete derived from the statistical procedure is a very important measure of the” 


yheric electric field, so is the spectrum derived in eo Se ee anothe: 
eeceert Brom nay point of view, whether an (oe! mospheric electri 


Se 
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The results are encouraging to understand the three dimentional structure of the 
atmospheric electric field, which may be regarded as a “turbulent field” consisting of 


various ranks of space charge concentration characterized by its own scale and its own 
density. 
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y—>a Transition in Fe203 with Pressure 


By Ikuo KUSHIRO 


Geological Institute, University of Tokyo 
(Read May 16, 1960; Received May 20, 1960) 


Abstract 


The temperature of transition of synthetic 7-Fe,0,; (maghemite) to a-Fe,O, 
(hematite) is lowered by increasing pressure. The transition line is represented 
by the equation, p=150—0.3T, where p is pressure in bar and J is temperature 
in °C. The transition of Ti-bearing 7-Fe,O, (titanomaghemite) needs higher pres- 
sure and temperature than that of pure 7-Fe,0,;. From the experimental results, 
it is suggested that maghemite is formed only near the surface of the earth, proba- 
bly at depths shallower than about 500m of the crust. 


Introduction 


The transition of cubic y-Fe,O; (maghemite) to rhombohedral @-Fe,0; (hematite) 

by heating has been examined by several investigators (Michel, 1937; Forrer, 1938; 
Chaudron and Michel, 1938; Pouillard, 1949; Bozorth, 1951 and Shibuya, 1951) and it is 
accepted that the transition occurs irreversibly at a certain temperature between 275 
and 800°C atlatm. The large difference in transition temperature has been interpreted 

as a result of presence of minor elements in solid solution. On the other hand, the 
effect of pressure on the transition temperature has not been known yet. In this paper, 
the results of the writer’s experiments on the y-a transition of Fe,O; in relation to pres- 
sure and temperature are presented. > 
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Table 1. Results of transition experiments 


eit aoe C y cre a Phases present 
690 415 23 hematite 
415 225 8 " 
BD 162 14.5 ' 
275 127 als}, 15) " 
, 160 351 12.5 " 
140 162 PAL aay " 
’ 105 130 1 maghemite 
70 313 215 hematite + maghemite 
70 250 24 maghemite 
Do 300 24 hematite -+ maghemite 
35 477 22 hematite 
35 236 10 maghemite 
20 300 10 " 
1 570 0.1 hematite 
al 530 24 hematite + maghemite 
1 490 15) maghemite 
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For comparison, the behavior of Ti-bearing 7-Fe,O; (titanomaghemite) was examin- 
ed. As the starting material, natural titanomaghemite from the Atumi dolerite, 
Yamagata Prefecture (Sp. No. 1003) was used. The composition determined by T. Katsura 
is as follows; FeO: 11.89, Fe,O3: 56.53, TiO, :. 22.63, MgO: 0.45, MnO: 2.29, V,O3: 0.73, 
Al,O3: 1.94 in wt. %. It is situated near the TiFeO, (ilmenite)-Fe,O, (hematite) join in 
the system TiO,-FeO-Fe,0O;. The lattice constant is 8.365+0.002 A. The experiment 
was done by the same method as that in the case of maghemite. Transition of the 
titanomaghemite to titanhematite did not occur even at higher pressure, but a decrease 
in lattice constant was observed. The results are listed in Table 2. The decrease in 
lattice constant, i.e. the decrease in unit cell volume, of titanomaghemite may be at- 
tributed to an existence of vacancies in the metallic position. This is supported by 
another experiment by which it was found that titanomagnetite from the Atumi dolerite 
(Sp. No. 1002) (FeO: 28:33 wt. %, Fe,O;: 44.18, TiO,: 18.10) having less vacancies was 
less compressed (8.410->8.402 A) than was No. 1003 at 2800 bars and 355°C. As shown 
in the table, the lattice constant ceases to decrease at higher pressure than 1200 bars, 
remaining almost at the same value of 8.344 A. This may indicate that the degree of 
compression attained to the maximum at about 1200 bars. It was also found that the 
titanomaghemite is stable up to higher temperatures than the stability temperature of 
pure maghemite. The specimen No. 1003 did not change at 675°C and even at 840°C 
at ordinary pressure. The stability of titanomaghemite is probably a function of the 
Ti content: higher pressure and temperature may be required for the transition of 
maghemite with higher Ti content. 


Table 2. Change in the lattice constant of titanomaghemite with pressure 


Lattice Apes 
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magnetite may indicate that the oxidation took place at temperatures higher than about 
900°C. 
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Propagation of Low Energy Cosmic Ray Particles Associated 
with Solar Flares 


By Kunitomo SAKURAI 


Geophysical Institute, Kyoto University 
(Read May 16, 1960; Received May 27, 1960) 


Abstract 

It has been observed that the low energy cosmic rays are produced associated 
with solar flares and propagate into interplanetary space. In this paper, we ex- 
amined the features of the propagation mechanism of these cosmic rays, taking 
into account the results deduced from the theory of geomagnetic effect on these 
low energy cosmic rays. In consequence, we can conclude that the features of 
their incidence to the earth must have the impact zone effects, G.M.T. effects, 
and seasonal variations. Then, we examined the relations between these low 
energy cosmic ray events and solar and terrestrial phenomena in view of the 


observed data and made clear that both of flares and type IV outbursts taking 


—— 


place on western side of the sun seem to correlate these cosmic ray events much 
better. Taking into account those results we tried to explain reasonably the 
propagation mechanism of these low energy cosmic rays into interplanetary space 
and in addition to estimate the properties of magnetic clouds in that space. 


1. Introduction 


During the period of I.G.Y., the observations of low energy cosmic rays associated 
with solar flares were performed and through these observations it has been made 
possible to deduce the electromagnetic state between the sun and the earth. Such low 
energy cosmic rays associated with solar flares was observed firstly at the time of the 
great solar cosmic ray event on February 23, 1956 and it was found that the low energy 
part of solar cosmic rays continued to decay with very long tail during a few ten hours — 
in the space around reams EE es 1956; van Allen and Winckler, 1956 ; Bailey, — 
ie 1959). eo during the period d of LGY., it was often observed that low 
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have some relations to occurrences of the solar radio emission type IV outbursts and 
the geomagnetic storms with sudden commencement (Hakura and Goh, 1959; Reid and 
Leinbach, 1959). And then, it has been made clear that the flares which occurred on 
the western side of solar surface have better correlations with the incidence to the 
earth of the low energy cosmic rays than with those on the eastern side (Reid and 
Leinbach, 1959). Moreover, Hakura, Takenoshita and Otsuki (1958) performed the detailed 
analyses about the events of low energy cosmic ray incidence observed on September 
13, 1957 and February 11, 1958 and made clear in detail the patterns of their incidence 
at the geomagnetic high latitude regions before the commencement and during the 
progression of geomagnetic storm. In order to investigate the propagation into inter- 
planetary space of the low energy cosmic rays from the sun, the manner of their in- 
cidence to the earth, and the electromagnetic processes occurring in the earth’s outer 
atmosphere in relation to the low energy cosmic ray events, such a work will be very 
important. 

In this paper, we shall examine the propagation mechanism of the low energy 
cosmic rays associated with solar flares making use of data published by Hakura and 
Goh (1959), Reid and Leinbach (1959) and other workers, which were obtained during 
the period from June, 1957 to May, 1959. In section 2, we shall at first investigate the 
features of the incidence to the geomagnetic high latitude regions of the low energy 
cosmic rays coming from the sun’s direction with the theory of geomagnetic effect on 
cosmic rays, and then show the results expected from the above calculations. Next, 
in section 3, we shall compare the results deduced from section 2 with the analyses by 
the authors quoted above, and examine the relationship between the incidence to the 
earth of such low energy cosmic rays and the solar and terrestrial phenomena associat- 


ed with them, and consider their propagation mechanism into interplanetary space. 


2. Motions of Low Energy Cosmic Rays in the Geomagnetic Dipole Field 


It is known that the solar cosmic ray events observed five times in the past show 
the large fraction of cosmic ray particles arrive directly at the earth’s neighbourhood 
from the sun and impinge upon the earth according to the expectation from theoretical 
calculations of the motions of cosmic rays in the geomagnetic dipole field, showing the 
so-called impact zone effects of their incidence (Firor, 1954; Liist and Simpson, 1957). 
The energies of these cosmic ray particles which produce such events range from about 
1 Bev to 30 Bev, but the energies of cosmic rays which we are going to examine here 
are much lower than the former, ranging from a few Mev to several hundred Mev. 
Thus, from the theory of geomagnetic effect of cosmic rays it is obvious that such 
low energy cosmic ray particles can only impinge upon the geomagnetic high latitude 
regions. 

In this section, referring to the results of calculations by Malmfors (1945), Dwight 
(1950), Schliiter (1951 ; 1958), Firor (1954), Stormer (1955), Jory (1956), and Liist (1956) 
and others and their methods of calculations, especially those of Stormer and Dwight, 
we shall examine the orbits of motion of cosmic ray particles with a few special rigi- 
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dities in the geomagnetic dipole field. The method of calculation is, roughly speaking, 


as follows; 

At first, we assume that a cosmic ray particle impinges vertically upon the geo- 
magnetic latitude 4. In Fig. 1, we shall show the directions of velocity vectors (Ao, 
¢) which a cosmic ray particle arriving at 4 will take before coming into the geo- 


magnetic field. 
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MAGNETIC RIGIDITY IN Gev 
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Fig. 2. Relation between the geomagnetic latitude of impact on the earth 


of particles arriving from the sun with magnetic rigidities of 30 Mev to 
3 Bev, and the geomagnetic latitude of source. 
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(2) At the higher latitude regions, since the geomagnetic axis does not coincide 
with the geographic one, it seems that a Greenwich time dependence of their incidence 
to the earth appears in addition to the effect (1). 

(3) The apparent position of the sun viewed from an observatory on the earth 
varies with the earth’s rotation and revolution. Thus, as the apparent position varies 
annually, the source direction of the low energy cosmic rays, which is now equivalent to A, 
the latitudinal component of velocity vectors, will vary in the same manner. Therefore, 
it may be expected that there is the seasonal variation for the incidence frequency of 
the low energy cosmic rays associated with solar flares. As an example for the above 
expectation, we shall show in Fig. 4 the case of American zone, the region near about 
geographic 70 degrees West meridian in the north hemisphere. Since | 4.| of hatched 
area is greater than 20 degrees, the low energy cosmic rays associated with solar flares 
will not be observable in the American zone if this zone is in that area. 
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particles are influenced very severely in their passage to the earth by the magnetic 
fields supposed to be existent in interplanetary space (see for example, Simpson, 1957; 
1958). Namely, it is regarded as very plausible that those interplanetary magnetic fields 
play an important role on the anisotropy degrees of solar cosmic ray incidence to the 
earth and make these solar cosmic rays to impinge upon the higher latitudes. There- 
fore, if we suppose that the propagation directions of such low energy cosmic ray 
particles as treated in this paper are disturbed considerably in the passage from the 
sun to the earth by the influence of those interplanetary magnetic fields, it may be 
expected that those low energy cosmic ray particles impinge isotropically upon the 
earth from the initial stage of their incidence. Although most of the low energy cosmic 
ray events continue for very long time, say several ten hours (Table I), it can not be 
expected that these long time durations of their incidence are produced by the con- 
tinuation of the solar production of low energy cosmic rays, because the associated 
flares continue at most during an hour. Consequently, there may be the mechanism 
which can trap these low energy cosmic rays for very long time in the solar neigh- 
bourhood or in the space between the sun and the earth and make them impinge upon 
the earth gradually. 


3. Associated Phenomena and the Propagation into Interplanetary Space 


If there is no magnetic field in interplanetary space, or if magnetic fields if any 
have configurations not affecting propagation of these low energy cosmic rays associated 
with solar flares, their incidence on the earth is to accord with the expectation obtained 
in section 2. So, we can expect that these low energy cosmic ray particles impinge 
predominantly upon the morning hemisphere of the earth. We may find out the above 
expectation in the work by Hakura et al. (1958), who investigated the low energy 
cosmic ray events on September 13, 1957. At present, since the investigations on the 
anisotropies of their incidence to the earth during a few hours after the beginning of 
their incidence have not been done in detail, we must say that our theoretical results 
concerning their incidence can not yet be confirmed. Therefore, in this section we 
shall not discuss in detail the features of their incidence to the earth, but examine 
the relationships among the productions of low energy cosmic rays at the sun, related 
solar phenomena and terrestrial ones, and thereafter, we consider the propagation 
mechanism of these cosmic rays into interplanetary space. 

A. Analyses of Observational Data 

Data used here were obtained in the period from June, 1957 to May, 1959 (Hakura, 
1959; Hakura and Goh, 1959; Reid and Leinbach, 1959). At first, we show eee 5 
the positions of solar flares connecting with the production of sae energy cosmic rays. 
In this figure, the type IV outbursts associated with flares producing these low energy 
cosmic rays and the ones not associated with these cosmic rays a0 also plotted together. 
We can see from this figure that the positions producing the incidence of low energy 
cosmic rays to the earth are mostly found on western side of the aii and thet sees 
flares are closely related with the occurrences of type IV outbursts and geomagnetic 
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Fig. 5. Distributions of the solar flares producing the low energy cosmic ray events, 
and of type IV outbursts, shown on the heliographic map. 
@: Flares associated with geomagnetic storms after low energy cosmic ray events. 
©: Flares not associated with geomagnetic storms after low energy cosmic ray 
events. : 
O: Type IV outbursts connected with low energy cosmic ray events. ; ? 
4: Type IV outbursts not associated with low energy cosmic ray boa 
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Table I, 


The Durations of the Impingement Upon the Earth of Low Energy Cosmic Rays 
(23 events are shown in this Table) 


Duration (hrs) Number of Events 
0O— 20 6 
21— 40 6 
41— 60 4 
61— 80 3 
81—100 4 


(a) 


AT 
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time of these cosmic rays becomes much longer than that of their direct transit. More- 
over, it is found from Fig. 7 that there is a tendency of the transit time becoming 
shorter and shorter as the positions of flares move more and more westwards on the 
sun. Thus, we can expect that these results obtained above give a very important clue 
to our consideration of the propagation mechanism of low energy cosmic rays through 
interplanetary space after their ejection from the sun. Finally, let us call attention to 
the continuation of their incidence during geomagnetic storms. 

B. Propagation Mechanism 

At present, it is well known that the cause of solar radio type IV outbursts is found 
in the relativistic electrons in solar outer corona (Boischot, 1957; Boischot and Denisse, 
1957). The existence of large number of these high energy electrons above several 
Mev leads us to deduce that many protons with energies of the same order are present 
there. Therefore, we can conclude that these protons with energies ranging from a 
few Mev to several hundred Mev are produced near the sun at about the same time 
as occurrences of type IV outbursts and then they begin to propagate into interplanetary 
space. Thus, when these high energy protons leave the sun, the plasma clouds moving 
in the solar outer corona may, too, be ejected from the sun together with these protons, 
and afterwards begin to propagate as so-called magnetic clouds. Although we can 
conclude that both of low energy cosmic rays and magnetic clouds are ejected from 
the sun together into interplanetary space and thereafter they arrive near the earth, it 
seems, aS we saw in subsection A, that there are some restrictions in the manner of 
the propagation of low energy cosmic rays. Therefore, taking into account these 
circumstances we must proceed to consider the mechanism in which low energy cosmic 
rays will already be affected very severely by the interplanetary magnetic fields before 
arriving near the earth since their incident directivities at the earth are very isotropical 
from the initial stage of their incidence (Reid and Leinbach, 1959). Both facts that 
the sun-earth transit time of them tends to being dependent on the heliographic positions 
of flares as can be seen from Fig. 7, and that as evident from Fig. 5 the heliographic 


positions of flares producing these low energy cosmic rays incline to western side of | 


the San, may give us the polaiing.< conclusion : namely, when arin eae rays 


‘ a Po —_— ae ee 


r , 


Propagation of Low Energy Cosmic Ray Particles Associated with Solar Flares 161 


periods of the solar activity maxima interplanetary magnetic fields may take such 
configurations as to disturbe violently the directivities of propagation of solar cosmic 
rays. So, referring to their investigation, we can suppose definitely the directivities of 
propagation of low energy cosmic rays treated in this paper through interplanetary 
space will be disturbed very violently by interplanetary magnetic fields. From the 
results obtained in this section and the consideration made by Towle and Lockwood 
(1959), we can deduce that there are magnetic irregularities affecting in the similar 
manner both of low energy cosmic rays with energies a few Mev- several hundred Mev 
and relatively high energy cosmic rays with energies of several Bev. It may be pos- 
sible to explain reasonably why Towle and Lockwood (1959) could not detect any solar 
cosmic rays with several Bev energies, if we presume that the interplanetary magnetic 
field configurations are variable according to solar activity cycles and therefore in the 
periods of solar activity maxima the incident directivities of these solar cosmic rays 
to the earth are purely random by the influence of magnetic fields mentioned above, 
and that the production rate, of these cosmic rays is much smaller than that of low 
energy cosmic rays. So, we may conclude with Towle et al. (1959) that since inter- 
planetary magnetic fields in the declining periods of solar activity seem to be transparent 
to the solar cosmic ray propagation, such results as Firor’s may be obtainable. Our 
explanation, however, is still unsatisfactory because the magnetic fields in interplanetary 
space deduced above indeed disturbe severely the incident directivities to the earth of 
low energy cosmic rays on the one hand, but on the other hand these fields have such 
properties as to maintain well their direction of propagation as a whole as shown in 
Fig. 5. Consequently, configurations of interplanetary magnetic fields must be able to 
explain these facts reasonably. If we refer as plausible to the magnetic fields in inter- 
planetary space presented by Cocconi et al. (1958) and Parker (1958a, b), it may be 
favorable to our reasonable explanation of the low energy cosmic ray events associated 
with solar flares. 

As evident in Table 1, most of the low energy cosmic ray events continue to 
impinge upon the earth during very long time and during geomagnetic storms, too 
(Freier, Ney and Winckler, 1959; Reid and Leinbach, 1959). Since it is very difficult 
to explain such long duration of their incidence on the basis of the solar flare origin 
alone, a mechanism has been proposed, in order to explain it reasonably, in which these 
low energy cosmic ray particles, as being trapped in magnetic clouds ejected from the 
sun, propagate into interplanetary space together with magnetic clouds (Gold, 1959). 
Thus it may be very plausible there is two manners in the propagation mechanism of 
low energy cosmic rays. Our explanation of the propagation mechanism of low energy 
cosmic rays will be as follows; particles arriving at the earth before the beginning of 
geomagnetic storms are trapped in magnetic clouds when leaving from the sun a 
+ diffuse out of them gradually, and most of the former propagate into 
interplanetary space along the radial magnetic fields of tongue type extended from the 
d those impinging upon the earth during geomagnetic storms arrive at the earth 
Moreover, the tail of their incidence remaining 
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after the end of geomagnetic storms may be produced by the leakage of particles 
trapped in magnetic irregularities near the sun. 

C. Characteristics of Magnetic Clouds 

It seems that magnetic clouds in which low energy cosmic rays are trapped will 
expand during the propagation into interplanetary space and their magnetic intensities 
will decrease gradually (Batchelor, 1955; Hirono, 1957). 

In this subsection, taking into acount that magnetic clouds are able to trap these 
cosmic rays for a long time, we shall examine, though very rough, the intensities of 
magnetic fields and their dimensions. Since we can suppose the cosmic ray particles 
impinging upon the earth in the initial stages of these events as having escaped from 
the trappings by magnetic clouds in a few 10 minutes or several hours, we can estimate, 
with those observations, the order of magnetic intensities or dimensions of young 
magnetic clouds near the sun. And moreover, taking into acount trappings of low 
energy cosmic rays for several 10 hours in later stages, we can estimate the lower 
limits on these quantities of magnetic clouds in interplanetary space. The results 
obtained with various magnetic intensities are shown in Tables II and III. The 


Table II. 
Dimensions of Magnetic Clouds near the Sun deduced from Low Energy Cosmic Ray Events 


Magnetic Field Dimensions of Magnetic Clouds (cm) 
Intensity : 
(Gauss) 0.5 hr* 1.0 hr Banos 2.0 DY es 2. 5 hr 
10-8 24x10% 35x10" 41x10" = 4.2x100 54x10" 
lOc PSOE 1010 d35<10! oe Oe 1.610" 


10s | oY : x ey 3. 5 x we 4. 1 x 10" : * 2x 10” 5.4x 10” 
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Lx=./Dr (cm), 


Where z + the time of escape above mentioned. Thus, we may obtain the dimensions 
of magnetic clouds, where in these estimations we assume E=10 Mev as representative. 

Although the dimensions of magnetic clouds in this subsection are much smaller 
than those considered usually at present, these values are near those estimated by 
Hirono (1957). 


4. Concluding Remarks 


In section 2, we examined the motions of low energy cosmic ray particles in the 
ee magnetic field of the earth and described some results expected from those 
examinations. Since the features of the incidence to the earth of those cosmic rays, 
however, have not yet been investigated in detail, we can not decide at present whether 
or not those particles show a tendency of impinging upon the earth in accordance with 
the expectations from theoretical calculations. We must, therefore, examine, in order to 
make clear the above mentioned problems, the manner of the incidence of these cosmic 
rays, taking into account observed data and referring to theoretical results. In section 
3, we discussed on both the solar and terrestrial phenomena occurring in connection 
with low energy cosmic rays associated with solar flares and showed that the arrival 
of these low energy cosmic rays at the earth has some relations with flares on western 
side of the sun, and that positions of type IV outbursts related with those events 
incline to western side. These facts seem to be favorable to the explanation that these 
low energy cosmic rays propagate into interplanetary space along the radial magnetic 
fields of tongue type extended from the sun and dragged somewhat by the solar 
rotation. It is estimated that in the propagations of these cosmic rays is included 
so-called indirect method propagation (Gold, 1959) of them trapped in magnetic clouds. 
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Abstract 


The daily percentage variations of the air-earth current were calculated every 
longitude of 90° over the world, using the percentage variation of the potential 
of the upper conducting layer (atmospheric total potential) deduced from Carnegie 
measurements over the oceans, and that of the columnar resistance obtained from 
the measurements by Sagalyn and Faucher. The results generally represent the 
world-wide distribution of the measured daily variation curves which were arranged 
by Israél. This fact generally indicates that the daily percentage variations of 
the columnar resistance at the most stations over the land area are roughly equal, 
and also of the same magnitude as the universal daily percentage variation of 
the atmospheric total potential. However, more detailed discussions by comparing 
the measured and the calculated results at the urban district of the meridian 135° 
E show that the daily variations of the air-earth current depend more on the 
local columnar resistance than the atmospheric total potential, the percentage 
variation of the former being twice as much magnitude as that of the latter. It 
is suggested that at the urban district the unclei concentration in the lowest layer 
of the exchange layer plays an important part in the atmospheric electric current 


system. 


1. Introduction 


The air-earth conduction current in the atmosphere has been considered to be of 
fundamental importance in the atmospheric electricity. The reason is as follows. First, 
the quasistatic state of the atmospheric electricity gives the constant value of the air- 
earth conduction current density with altitude from the upper conducting layer to the 
earth; thus the time variation of the air-earth current has been expected to be an 
important clue to indicate the time variation of the potential of the upper conducting 
layer with respect to the earth. Second, the air-earth current depends on the atmos- 
pheric total potential and the columnar resistance and they will be both less subjected 
to variations than the potential gradient and the conductivity ; so the air-earth current 
will be less subjected to variations with time of day. Thus, the air-earth current has 
also been expected to be less influenced by local effects; local changes of conductivity 
can alter the air-earth current only to the extent to which they alter the tolal co- 
lumnar resistance. These should be experimentally confirmed. 

Although the statistical method has shown the characteristic definite types of daily 
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variation of the potential gradient, the air-earth current does not seem to tend toward 
definite types of variation but give rather divergent character. Miihleisen (1956) pointed 
out that the air-earth current is greater in the daytime than at night, in spite of the 
increasing of the columnar resistance in the daytime. He has concluded that the posi- 
tive space charges which enlarge the primary field in the lowest air layers must be 
considered. Though the space charge effects must be certainly taken into consider- 
ations, they do not cause the divergent character of variation of the air-earth current 
only by themselves. 

The daily variation of the potential of the upper conducting layer depends on the 
Greenwich Mean Time and that of the columnar resistance depends on the local time. 
As the air-earth current is related to those two factors, the daily variation of the air- 
earth current should depend on the longitude. Therefore, the fundamental types of 
variation to be in any longitude of the world should be considered. And also it should 
be examined that which the daily variation of the air-earth current at every location 
of the world depends on more the atmospheric total potential or the columnar re- 


sistance. 


2. Longitude Distribution of the Types of Diurnal Variation of the Air-Earth Current 


The air-earth current can be measured by two different methods. The actual air- 
earth current can be measured directly, and also it is possible to give the air-earth 
conduction current indirectly by measuring the local potential gradient and the local 
conductivity independently. The air-earth conduction current density i is related to the 
potential gradient E and the conductivity A by Mayan 

t=A Bs ’ ee 


On the other hand, the air-earth conduction current density is also given by 
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by the percentage variations of E and 4, or V and R without considering their abso- 
lute values. 
Israél (1954) has arranged the diurnal curves of the air-earth current at 10 stations 
in the world according to the geographic latitude, which are reproduced?int Figs 
Between all of them he could not find out any common feature depending on either 
the universal time or the local time. Thus the systematical explanation of such diurnal 
curves has not yet been given. 


0 6 IZ ol8a 324 


LOCAL TIME LOCAL TIME 

Fig. 1. Mean diurnal variations of air-earth current density in percentage (obtained 
during the period of one year at least) at the stations 

a) Fairbanks/Alaska (Gish and Sherman, 1940) 

b) Uppsala (Benndorf and Hess, 1928) 

c) Potsdam (Kahler, 1909) 

d) Kew (Scrase, 1933) 

e) Chambon-la-Forét (Institute de Physique du Globe, Paris) 

f) Buchau (Israél, 1954) 

g) Tucson (Wait and Torreson, 1941) 

h) Huankayo (Peru) (Torreson and Wait, 1948) 

i) Bandoeng/Java (Leckie, 1938) 

k) Watheroo/Westaustralien (Wait and Torreson, 1941). 
The stations were arranged according to the geographic latitude. 
abscissas give 0% GMT. (after Israél, 1954) 


The arrows on the 


In the equation (4), if the percentage variation of FR is negligibly small compared 


with that of V, the percentage variation of 7 is the same as that of V and must depend 


on the Greenwich Mean Time. While, if the percentage variation of V is negligibly 


small compared with that of R, the percentage variation of 7 undergoes the change in 
the opposite sense to R and must depend on the local time. However the world-wide 
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distribution of the diurnal curves of the air-earth current given by Israél indicates that 
neither the percentage variation of the atmospheric total potential nor that of the co- 
lumnar resistance are negligible and the two quantities are both quite effective on the 
variation of the air-earth current; the author sees the coupling effect of the two quanti- 
ties in the Israél’s empirical results. 

The potential of the upper conducting layer cannot be directly measured, but it is 
possible to estimate the potential from the measured potential gradient distribution 
with altitude by integrating from the earth to the upper conducting layer. This is 
given by 


v=—f Bah, (5) 


where H is the height of the presumed upper conducting layer. Clark (1958) estimated 
the mean potential of 290 kV from the measurements between 0.015 km and 6km above 
locations from Greenland to Key West by an airplane. 

From the equations (1) and (2), the potential gradient is shown by 


=1/ (6) 


By differentiating (6), the time variation of the potential gradient is given by 
IdE Ji-dV isl dd, idk 
—s (7) 


E dt V dt i dt. Ride 


If 2 and R are unaffected Py time of day, then E shows the variation of V. These 
conditions are achieved over the ocean. 
In analysing the results obtained during 
the cruise of the Carnegie, Mauchly 
(1923) showed the variation of potential 
gradient over the ocean to be singly 
wpenodics a a BiesaaTaDS at about at 


of the local ti 
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and Faucher (1956). 


are given in Fig. 3. 


sistance in any other locations over the l 


The daily percentage variation of the average columnar resistance 


It is not proved whether the daily changes in the columnar re- 


In percentage as that in New Hampshire. 


and area are of the same order of magnitude 


Using the percentage variation of the atmospheric total potential in Fig. 2 and 


that of the col r resistance in Fig. ‘ iati 
ummnar resistance in Fig. 3, the diurnal percentage variations of the air- 


earth current can be calculated by the equation (4). The results are obtained every 


LOCAL TIME 


Fig. 4. Diurnal variations of air-earth current 
calculated on every longitude of 90°. 


give 0” GMT. 


Table 1. Relations between the calculated and the measured air-earth currents 
(Circles indicate good agreement and cross marks indicate poor agreement 
in the general features of the period and the amplitude.) 


The arrows 


longitude of 90° and shown in Fig. 
4, in which the arrows give 0* GMT. 
It is interesting to see how repre- 
sentative the characteristic features 
in these diurnal curves of variation 
over the world are of those in the 
measured ones given in Fig. 1. The 
most characteristic feature in Fig. 4 
is that there are two types of diurnal 
variation in the air-earth current. 
One is the double periodic variation 
with one maximum before and after 
noon respectively—curve A. The 
other type of variation is a single 
periodic one with a maximum in the 
morning and a minimum in the after 
noon—curves 8, C and D. The 
double periodic type of variation has 
smaller amplitude in daily change 
than the single periodic one. 

In comparison between Figs. 1 
and 4, the relations between the 
measured and the calculated vari- 
ation curves of the air-earth current 
are shown in Table 1, in which the 


Calculated A By thee D 
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circles indicate good agreement between them in the general features of the period 
and the amplitude, while the cross marks indicate poor agreement. As the measured 
air-earth current may contain the effects of other generators than the thunderstorm 
generator and the columnar resistance used in the calculations may require some 
corrections at the respective stations, it is difficult to expect more excellent agreement 
between the two cases. From the general agreement of the calculated and measured 
air-earth currents, however, it can be pointed out that the diurnal percentage variations 
of the columnar resistance in any locations over the land area are roughly of the same 
order of magnitude as that in New Hampshire, and therefore the columnar resistance 
as well as the atmospheric total potential has an important role in the atmospheric 


electric current system over the land area. 


3. An Example at the Urban District of the Meridian 135°# 


In order to give more detailed discussions, the diurnal percentage variation of the 
air-earth current is calculated at the meridian 135° (the standard meridian of Japan), 
using the above percentage variations of the atmospheric total potential and the co- 
lumnar resistance in Figs. 2 and 3. The result is shown in Fig. 5, in which the time 
is referred to the mean local time of the meridian 135°L. 

At the Geophysical Institute of Kyoto University, the potential gradient, polar 
conductivities and the wind velocity have been recorded on a recording paper. The 
Institute is located in the north-east district in Kyoto City which has the population 
of about one million and no big industrial area. Throughout the year and especially 
in the cold season, it is excelled in westerly wind over Kyoto. The generating volt- 
meter type field meter and the Gerdien type conductivity meters with the vibrating- 
reed electrometers are set atop the building 15m high to avoid disturbances near the 
ground. Using the values of the measured potential gradient and conductivities, hourly 
values of the air-earth current are calculated on every day. 

Although the potential gradient and the conductivity are largely subjected to 
local effects, they appear usually as a mirror image of the time variation of each 
other; so their product, that is, current density will be generally expected to be > more — 
Ne ss amversal eae the rebar current res obtained are e 
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LOCAL TIME 
Fig. 5. Diurnal variation of air-earth current 


calculated in the longitude of 135°E. The arrow 
gives 0% GMT. 


0 4 8 12 16 20 24" 


LOCAL TIME 
Fig. 6. Average diurnal variation of air-earth 
current obtained indirectly in Kyoto. The arrow 
gives 0% GMT. 


on every month are calculated. In 
Fig. 6, the average curve of variation 
of 4 months from September to De- 
cember, 1959 is shown. 

The most remarkable features in 
Fig. 6 are that the amplitude of the 
variation is about twice as that of 
the fundamental type of variation in 
Fig. 5, and the time of appearance 
of the minimum displaces toward 
later. Such differences can be sum- 
marized in terms of the columnar 
resistance in the atmospheric electric 
current system over the land area. 
Using the percentage variation of the 
atmospheric total potential in Fig. 2, 
the diurnal percentage variation of 
the columnar resistance can be calcu- 
lated. The result is shown in Fig. 7, 
which represents the characteristic 
features of the columnar resistance 
at the urban district. 

The result is compared with the 
columnar resistance in Fig. 3. In the 
first place, the amplitude of the vari- 
ation is about twice as that in Fig. 
3. This shows the daily variation of 
the columnar resistance to be about 
twice as much magnitude as that of 
the atmospheric total potential. Sec- 
ond, the time of occurrence of the 


minimum in the morning is 4” and the substantial increasing begins at about 6”, while 
the minimum appears at 7” in Fig. 3. This difference will depend on the degree of 
the urban conditions over the land area. In urban districts, the increasing of the co- 
lumnar resistance is caused by the accumulation of the unciei in the lowest layer of 
the exchange layer, while in rural districts, it is the result of conveyance of the nuclei 
from the urban districts. It results in the difference of the time of the occurrence of 
the increasing of the columnar resistance. Third, the time of occurrence of the maxi- 
mum in the afternoon is 18”, while it appears at 15" in Fig. 3. This is related to the 


mechanism of the daily change in the atmospheric electric field, which will be described 


in a later paper. ; ; 
Early in the morning before the sunrise when atmospheric state is the most stable, 
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in urban districts, the conductivity at 


130 the earth’s surface begins to decrease 
120 by the production of man-made nuclei 
110 which accumulate in the lowest layer 
100 of the exchange layer. It results in the 
90 increasing of the columnar resistance. 
80 In the daytime when the convection be- 
70 comes stronger and the wind begins to 
60 blow, the turbulent mixing occurs and 

ns destroys the nuclei accumulation in the 

ie) 4 8 12 16 20 24 


lowest layer. The nuclei are brought 


eS up into the upper layers in the ex- 
Fig. 7. Average diurnal variation of columnar 
resistance deduced from Figs. 2 and 6. Dashed 
line represents the case without horizontal diffu- 


sion of nuclei. The arrow gives 0% GMT. results slackens the increasing of the 


change layer, and at the same time 
they are diffused horizontally. As the 


total resistance of an air column of 
unit area in horizontal cross section. The effect of this mechanism can be seen from 
9" to 17” in Fig. 7. The columnar resistance presumed to be unaffected by the hori- 
zontal diffusion of nuclei is represented by the dashed line in the figure. In comparison 
between Figs. 3 and 7, the maximum at 15” in Fig. 3 falls in the interval of this dashed 
line in Fig. 7; the maximum of the columnar resistance in rural districts depends on 
the horizontal conveyance of the nuclei produced at the sources in urban districts. 
The distribution of daily variation of the columnar resistance with distance from the 
sources of nuclei should be investigated. 


4. Concluding Remarks — 


According to Kraakevik (1958), the air-earth conduction current density in the ex- 
change layer is larger than that above the exchange layer on the average some 30%, 
due to upward convection of positive space charge. This is the effect of the local or 
the uence generator. prot this generator Miihleisen (1958) has investigated in 
details, aa it i difficult to” to separate the ie of ee local ees Bei: 
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